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Abstract. A strong light pulse propagating in a nonlinear medium causes an effective change in the local 
refractive index. With a suitable tuning of the pulse velocity, the leading and trailing edge of the pulse 
were predicted to behave as analogue black and white horizons in the limit of a dispersionless medium. In 
this paper, we study a more realistic situation where the frequency dispersion of the medium is fully taken 
into account. As soon as positive frequency modes with negative norm are present in the comoving frame, 
spontaneous emission of quantum vacuum radiation is expected to arise independently of the presence 
of horizons. We finally investigate the kinematic constraints put on the emission and we show that the 
optimal directions to observe Hawking-like emission form a narrow angle with the direction of propagation 
of the pulse. 

PACS. 04.62.+V Quantum fields in curved spacetime - 42.65.-k Nonlinear optics 

1 Introduction 

Hawking radiation [1,2] is the quantum production of particles from vacuum fluctuations due to the presence of a 
black hole horizon in a curved stationary geometry. However, this prediction does not reside on the peculiar dynamical 
features of a general-relativistic spacetime, but only on the kinematic properties of a quantum field living in a curved 
spacetime. Following the pioneering work [3], analogous quantum vacuum emission phenomena have been predicted to 
occur in other physical systems, provided that a quantum field propagates on an effectively curved spacetime endowed 
with a horizon (see [4] for a complete review on analogue models of gravity). 

In particular, optical systems provide promising experimental analogues of curved spacetimes [5,6], because they 
are quite simple to handle and the technology to perform the required measurements is well developed. In fact, the first 
claim of analogue Hawking radiation in a laboratory has been recently reported in [7, 8] , where the analogue horizons are 
created by sending a strong laser pulse through a nonlinear dielectric medium with an appropriate velocity. The claimed 
observation of optical Hawking radiation was obtained by looking perpendicularly to the direction of propagation of 
the pulse, so to eliminate other spurious effect. Unfortunately, this result is still considered as controversial by some 
authors, who recently raised a few issues [9,10]. Nonetheless, no convincing alternative explanation of the experimental 
observations has been found yet. It is worth mentioning some attempts to interpret them as a dynamical Casimir 
emission [11] and/or the optical analogue of cosmological particle creation [12]. 

In this paper, we extend the existing models for analogue Hawking radiation produced by a strong pulse propagating 
through a nonlinear optical medium, by considering the kinematic effect of the frequency dispersion on the emission 
process. This allows to put constraints on the observable features of the emission and suggests the most suitable 
experimental conditions for its observation. 

Focusing our attention on a fused silica glass as the nonlinear medium, we first study a one-dimensional system. In a 
nondispersive stationary medium, the necessary and sufficient condition to trigger particle creation a la Hawking would 
be the presence of analogue horizons and the spectrum of the spontaneously emitted particles would be Planckian 
at a temperature proportional to the surface gravity of the horizon [13]. In dispersive media such as fused silica, the 
situation is kinematically more complicate, but quantum vacuum emission remains still possible: while horizons are no 
longer necessary for the emission, a sufficient condition is the presence of positive-norm modes with negative frequency 
in the reference frame comoving with the pulse. However, in this case, there is no a priori reason for the spectrum to 
be thermal. 

In the second part of the paper, we study the propagation of photons in a three-dimensional system with axial 
symmetry. In particular, we investigate how the kinematical properties of the emission in the comoving frame transfer 
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into the observable radiation in the laboratory and we determine which components of the quantum vacuum radiation 
are able to actually exit the glass reaching the detector. 



2 Mode analysis in one-dimensional media 

We consider a laser pulse with frequency fl p , propagating with velocity V in a parallelepiped of glass. If the pulse is 
strong enough, the refractive index experienced by a probe field increases [14] proportionally to the laser intensity /: 

n(I) = no + Sn(I) = no + n-2 1, (1) 

where n 2 «3x lCT 16 cm 2 /W. For a realistic value of I « 3 X 10 12 W/cm 2 , 6n » 0.001 [7,8]. 

Consequently, the velocity of propagation v- m of the probe field inside the pulse is smaller than outside (v ou t)- When 
neglecting the effects of dispersion and if the parameters are arranged such that v- ln < V < v out , photons propagating 
inside the pulse cannot cross the leading edge, but they are all dragged towards the trailing one. The leading (trailing) 
edge is then the analogue of a black (white) horizon, as seen by an observer outside the pulse [5]. Since Hawking's 
particle production depends only on the dynamics of the probe field in the analogue geometry, the presence of an 
analogue black horizon is a sufficient condition to trigger the Hawking process. 

In the presence of dispersion, such a simple picture is more complicated. In dispersive media, indeed, quantum 
vacuum emission is possible even without analogue horizons, as shown in [15] in an analogue model based on a Bose- 
Einstein condensate. In this section, we shall explore both horizon and horizonless configurations, where photons can 
be spontaneously emitted in a one-dimensional model of glass perturbed by a strong laser pulse. For definiteness we 
shall work with fused silica. The dispersion relation in such a medium relates the wave number K, measured inside 
the glass, to the frequency J7: 

c 2 K 2 =nl(n)Q 2 , (2) 
where the refractive index is given by the Sellmeier relation [16] 

(n/n,) 



with 

Ci = 0.8974794, hf2 1 = 0.125285 eV, 

C 2 = 0.4079426, Hf2 2 = 10.6661 eV, (4) 

C 3 = 0.6961663, hf2 3 = 18.1252 eV. 

In fig. 1, left panel, we plot all the four branches of the dispersion relation, but, in what follows, the analysis is 
restricted to the optical branch (red line), between the upper edge of the lowest gap at H(2 gap — 0.149636 eV and the 
next pole hfl-i. We give below reasons for this choice. 

To proceed, it is convenient to write the dispersion relation in the frame comoving with the pulse 



c 2 



kl+^ikx + Vu/c 2 )' = 1 2 (u; + Vk x fn( 1 [u J + Vk x },I)\ (5) 



by applying to eq. (2) a Lorentz boost at velocity V 

n = j(v)[u) + vk x ], 

K x = 1 (V)[k x + Vu/c\ (6) 

K — h 

J\y — fry. 

Hereafter, lowercase (uppercase, respectively) will denote quantities measured in the comoving (glass rest, respectively) 
frame. In fig. 1, right panel, the optical branch (red) and the lowest branch (blue) of the dispersion relation are drawn 
in the glass rest frame (K x ,fi). Positive- and negative-frequency branches are plotted with thick and thin lines, 
respectively. The coordinate system is then transformed under a boost with V/c — 0.66. This value of the velocity is 
chosen to realize a configuration with an analogue black hole horizon as it will be shown in the following section. In 
practice, it can be obtained as the group velocity of a pump laser pulse with a wavelength of approximately 400 nm. 
In what follows, we work with frequencies much lower than i?2, such that branches above the optical one cannot be 
excited. We shall therefore neglect the contribution of those branches, which are not shown in the right panel of fig. 1. 

In this paper, we discuss only the kinematic properties of light propagation in the considered system. The analysis 
of the spectrum of the vacuum emission is postponed to further investigations [17]. 
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Fig. 1. Left panel: Dispersion relation of fused silica (2) and zoom (inset) on the small frequency region. The red line is 
the optical branch. The blue line, hardly visible in the main plot, is the lowermost branch. Right panel: Optical (red) and 
lowest (blue) branches of the dispersion relation (2) in the glass rest-frame. The new axes {k x ,uj) are obtained by applying a 
Lorentz transformation with V = 0.66c to the axes (K x , fi). Thick and thin lines are positive- and negative-frequency branches, 
respectively. 



2.1 Horizon configurations 



We start by considering the configuration shown in fig. 2. To identify the solutions of the dispersion relation for a given 
comoving frequency cj, in the figure we plot the lowest (blue) and optical (red) branch of the dispersion relation in 
the frame of the pulse for the same pulse velocity V as in fig. 1, right panel. Positive (negative) laboratory frequency 
branches are represented by thick (thin) curves. They correspond, respectively, to positive- and negative-norm modes 1 . 
For illustrative purposes, we choose an extremely high (possibly not so realistic) intensity / « 3 x 10 14 W/cm 2 , yielding 
a strong perturbation of the refractive index Sn = 0.1. The dispersion relation inside the pulse (where Sn = 0.1) and 
outside the pulse (where Sn = 0) are plotted in the left and right panel, respectively. The top panel reports the 
positions where the dispersion relation is computed, with respect to the leading part of the pulse. We describe the 
structure of the scattering modes only around this leading edge. The discussion of the trailing edge follows the same 
lines. 

From fig. 1, right panel, one sees that the lowest branch (blue line) saturates to fl\ at small wave numbers. Except 
for very low comoving frequencies w (see fig. 2), excitations on this lowest branch have a wave number large enough 
that the rest-frame frequency f2 is equal to fl\ , independently of the momentum K x . As a consequence, these particles 
cannot propagate, because they have zero group velocity in the glass rest frame. Eventually they are completely 
absorbed. Thus, we shall consider only the optical branch. The case in which w is so small that the lowest branch 
cannot be neglected requires a different analysis, beyond the scope of this paper. 

We now introduce the concept of frequency-dependent horizon. Assuming that the pulse is moving in the positive 
x direction (V > 0), such a horizon is present for a frequency uj when inside the pulse there are only negative 
group velocity (measured in the comoving frame) modes, while outside there are both negative and positive group 
velocity modes. In this case, light can propagate only left-ward inside the pulse, both left- and right-ward outside it. 
Equivalently, there exists a point in between the internal and the external regions such that the speed of light in the 
comoving frame vanishes at this frequency, 

£=°- « 7 » 

In fig. 2, the dispersion relation is graphically solved on the optical branch (red curves) for a generic frequency 
(solid horizontal line). The arrows indicate the direction of propagation of the corresponding modes, which are labeled 
in the left and right panel by x (0) if they have negative (positive) laboratory frequency f2, or, equivalently, negative 
(positive) norm. The subscript in (out) denotes incoming (outgoing) modes, having group velocity directed toward 
(away from) the horizon. Note how the nontrivial shape of the dispersion relation (2) is responsible for new interesting 
features with respect to the standard models of subluminal dispersion usually considered in the literature [4]. 

For the chosen value of the comoving frequency lu, eq. (7) is satisfied at a point where Sn w 0.017. For this particular 
value of Sn, the dispersion relation (central panel) is tangent to the line of constant w (two coincident solutions denoted 



1 A rigorous proof of this fact will be provided in [17], by using the norm naturally induced by the Hamiltonian structure. 
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Fig. 2. Optical branch (red) and lowest branch (blue) of the dispersion relation (5) in fused silica for vanishing transverse 
momentum k y , in the frame comoving with a pulse moving with velocity V = 0.66c. Positive (negative) laboratory frequency 
branches are represented by thick (thin) curves. The dispersion is plotted for three values of the perturbation of the refractive 
index: Sn = 0.1 (left panel), Sn = 0.017 (central panel) and Sn = (right panel). The black horizontal line represents a generic 
frequency for which there are 4 real solutions in the external region (right panel), 2 real solutions in the internal one (left 
panel), and a point where the solutions passes from 4 to 2 (central panel). In the right panel, the dashed horizontal lines are 
the maximum and the minimum values of the frequency for which there are 4 solutions in this external region. The top panel 
reports the positions at which the dispersion relation is computed, with respect to the leading part of the pulse, around the 
analogue black horizon. In the left and right panels, the arrows indicate the direction of propagations of modes corresponding to 
solutions on the optical branch. Modes with negative (positive) norm or, equivalently, negative (positive) laboratory frequency 
O, are labeled by x (4>)- Subscript in (out) refers to incoming (outgoing) modes. 



by an empty dot). For outer points (right panel, Sn < 0.017), the dispersion relation has 4 solutions. Three of them 
correspond to incoming modes (Xin, <fc^, and since they propagate from x = +oo toward the horizon, having 
negative group velocity. The fourth one (denoted by an empty dot) corresponds instead to an outgoing mode ((f>^ ut ), 
with positive group velocity. For inner points (left panel, Sn > 0.017) the dispersion relation has only 2 real solutions, 
both corresponding to outgoing modes (xout and 4>o Ut ), which, having negative group velocity, propagate from the 
horizon toward a; = -oo. As a result, according to the definition given above eq. (7), within the frequency range 
u min < ui < ui max (where hcj m i n ps 0.034 eV and fiw max ps 0.123 eV are represented in the right panel with dotted lines), 
the leading edge behaves as a black hole horizon. Analogously, the trailing edge behaves as a white hole horizon. In a 
more realistic situation, for typical experimental parameters [7], one has a much smaller Sn. As a result obtaining a 
configuration with horizon requires a very accurate tuning of the pulse velocity, down to the fourth digit. For instance, 
with Sn ps 0.001, a suitable value of the pulse velocity is v ps 0.6838c. Furthermore, the horizon turns out to be present 
only in a very narrow frequency window, between Hoj m i n ps 0.01379 eV and hui max ps 0.01396 eV, for the considered 
parameters. 

Summarizing, for all frequencies w m in < w < w max , there are 3 incoming (xm, <j>^ and 0^) and 3 outgoing {Xouti 
4>o U f an d <^out) modes. Consequently, the scattering (conversion of incoming modes into outgoing ones) at the analogue 
black horizon is described by a 3 x 3 matrix [18]. For instance, if one sent a negative-norm wave Xin into the system, 
it would be decomposed in a superposition of the three outgoing modes (see fig. 3), 

Xin — > aXout + 50out + Wouf (8) 

The three coefficients of this decomposition form a row of the scattering matrix. Analogously, the remaining two rows 
are obtained by decomposing the other two incoming modes. Since the scattering matrix mixes modes with positive 
and negative norm/laboratory frequency, the initial vacuum (i.e., the quantum state with no incoming particles) is not 
equivalent to the final vacuum. As usual, this implies that quanta on outgoing modes can be spontaneously emitted 
starting from initial vacuum fluctuations. In particular, it is easy to show that the occupation number of the positive- 
norm outgoing mode <^ ut is given by |/3| 2 . In a different analogue model based on flowing atomic condensates [19], 
scattering is also described by a 3 x 3 matrix and the spectrum of spontaneously emitted particles, given by |/3| , 
follows a Planckian distribution, as for standard Hawking radiation, but only up to a certain frequency w max [20,21, 
22]. Shall we observe Hawking-like radiation also in the present system? Unfortunately, the scattering is now described 
by a 3 x 3 matrix only for uj > uj m i n . In fact, the two modes (j>^ ut and (j)f n disappear for u < w m i n , thus the scattering 
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Fig. 3. Decomposition of the incoming mode Xin into outgoing modes. This mode travels from x = +00 toward the horizon 
where it is scattered in the three outgoing modes Xout, </>out> an d 4>out- Positive (negative) norm modes are represented by solid 
(dashed) lines. 

matrix becomes 2x2 and (3 is not defined. However, since a negative frequency incoming mode (xin) is still present, 
quantum vacuum emission is expected to appear also in this case, albeit with different spectral properties. 

2.2 Horizonless configurations: fast pulse 

As a second case, we consider a pulse propagating with a large velocity V, such that the probe group velocity v is smaller 
than V for any frequency u). fn this situation, no positive group velocity modes propagate cither inside or outside 
the pulse, i.e. condition (7) is not satisfied. This configuration was experimentally realized in [7,8] and theoretically 
investigated in [23], where the nonlinearity of the dispersion relation is neglected in the quantization of the probe field. 
Dispersive effects have also been taken into account in a recent analysis [24] . In fig. 4, we plot the dispersion relation 
inside the pulse (left panel) and outside (right panel) for this experimental set up (V — 0.69c, Sn — 0.001). Thick 
(thin) lines denote positive (negative) norm/laboratory frequency branches. With the same convention introduced in 
the above section, modes with negative (positive) norm are named x (</>)• The dispersion relation has now only two 
real solutions on both sides. All of them correspond to leftgoing modes, two of them correspond to incoming modes 
Xin and O ut> and two of them to outgoing modes Xout and </>out- Since the dispersion relation has now only two real 
solutions, the scattering process is described by a 2 x 2 matrix. An incoming negative laboratory frequency mode Xin 
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Fig. 4. Positive (thick curves) and negative (thin curves) norm/laboratory frequency branches of the dispersion relation (5) 
in fused silica for vanishing transverse momentum k y , in the frame comoving with a pulse moving with velocity V ~ 0.69c. 
Sn = 0.001 (left panel) and 5n = (right panel) as in the experiment of [7,8]. In both the internal (left) and external (right) 
regions only two real solutions are present. Arrows indicate the direction of propagations of modes. Modes with negative 
(positive) norm or, equivalently, negative (positive) laboratory frequency Q, are labeled by \ ((f)). Subscript in (out) refers to 
incoming (outgoing) modes. 
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Fig. 5. Positive (thick curves) and negative (thin curves) norm/laboratory frequency branches of the dispersion relation (5) 
in fused silica for vanishing transverse momentum k y , in the frame comoving with a pulse moving with velocity V ~ 0.66c. 
Sn = 0.015 (left panel) and 5n — (right panel). There exists a threshold frequency cj t (left panel) such that, for uj t < oj < U) m&K , 
the number of propagating modes is 4 in the external region (right) and 2 in the internal one (left), as in the system of fig. 2. 
However, for uj < uj < uj t (horizontal solid line), the number of modes is 4 in both regions, similar to what was found in [15,25] 
for superluminal (supersonic, for analogue BEC systems) systems with superluminal (supersonic) dispersion relation. Finally, 
for uj < o; m i n oru> CJ ma x, there are two real solutions of the dispersion relation. Arrows indicate the direction of propagations 
of modes. Modes with negative (positive) norm or, equivalently, negative (positive) laboratory frequency O, are labeled by \ 
(<f>). Subscript in (out) refers to incoming (outgoing) modes. 



is still present and its decomposition in the two outgoing modes is 

Xin > &Xout + B(j) out . (9) 

The disappearance of the horizon and, consequently, of the outgoing mode </)^ ut removes the Hawking-like channel of 
particle production. Nevertheless, particles can be spontaneously emitted because the scattering matrix still mixes 
positive- and negative- norm modes. For instance, the occupation number of the outgoing mode </) out is given by \B\ 2 . 



2.3 Horizonless configurations: slow pulse 

As the last case we consider a pulse which is slow enough that, for some frequencies, the group velocity v of the probe 
field is larger than the pulse velocity V, both inside and outside the pulse. This implies that no horizon is present for 
those frequencies. Nevertherless, for all values of uj, negative-frequency modes are present, so that quantum vacuum 
emission is expected. This is possible because the glass rest frame frequency Q goes to a finite value Q2 for k — > 00. 

In this configuration, the dispersion relation (fig. 5) has now a maximum both in the external (dn ~ 0, right panel) 
and in the internal region (8n = 0.015, left panel), for a pulse velocity V = 0.66c. Consequently, there is a threshold 
frequency Wt (left panel), such that, for uj t < uj < w max , the system is described by a 3 x 3 scattering process as in 
configurations with horizons. For uj < uj m - ln and for uj > uj ma ,x, there are instead only 2 real solutions in both regions, as 
in the above discussed case. Finally, for w m ; n < uj < uJt (horizontal solid line), there are 4 real solutions in both regions. 
They correspond to 4 incoming modes, with group velocity directed toward the horizon (<j>f n , propagating rightward in 
the left region; Xin, <?^j and (j)f n , propagating leftward in the right region), and 4 outgoing modes, with group velocity 
directed away from the horizon (xout, 0out; an d <^out: propagating leftward in the left region; 0?, propagating rightward 
in the right region). Two (x; n and Xout) of these 8 modes have negative norm/laboratory frequency. As a consequence, 
the scattering is now described by a 4 x 4 matrix, and the decomposition of the incoming negative-frequency mode is 

Xin — ► aXout + £-0" ut + B + 0+ u t + /^ouf (10) 

From this expression one can read the occupation numbers |-B - | 2 , |-B + | 2 , and |/3| 2 of the positive-norm outgoing modes 
(f)~ ut , 4>o Ut , and 0^ ut , computed on the initial vacuum of the incoming modes. A similar situation was studied in [15, 
25] in a Bose-Einstein condensate, where the dispersion relation is supersonic. In that case it has been shown that 
radiation is spontaneously emitted, but the spectrum is not Planckian. We expect that those results can be directly 
extended to the present situation since, for a given comoving frequency ui, a system with subluminal dispersion relation 
where v > V is equivalent to a system with superluminal dispersion relation where v < V [26] . 
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Fig. 6. Loci of constant emission frequency u) (solid lines), equally spaced in log(cj), in the comoving frame (k x ,k y ), with 
V — 0.66c, inside the glass. The dashed lines are uj = uj m i u (blue), and ui = oj ma x (red). The cyan curve refers to a frequency 
uj < ojmin, from blue to red to o; m ; n < ui < oj ma x, and from red to yellow to uj > oj max . In the right panel the low momentum 
region is zoomed. A mode with positive (negative) norm/laboratory frequency does not undergo total reflection and can cross 
the interface between glass and air, if the corresponding wavevector (k x ,k y ) is contained in the green (blue) area. 

3 Propagation in three-dimensional systems 

In the above section we described the configurations where quantum vacuum emission may be present in an optical 
system. For the sake of simplicity we restricted our analysis to a one-dimensional model, with light propagating on 
the axis parallel to the pulse velocity. Focusing on the configuration with analogue horizons (sec. 2.1), we now relate 
the features of the emission in the comoving frame to what can be detected in the laboratory, outside the glass. 
We therefore consider the propagation of emitted photons in a three-dimensional system with axial symmetry, thus 
effectively two-dimensional. Moreover, we assume that the direction x of propagation of the pulse is parallel to the 
interface between glass and air, through which observations are made, as in [7]. The y-axis is then orthogonal to this 
interface. For the sake of simplicity, the air refractive index is taken equal to 1. Frequencies and wave numbers will be 
denoted by the subscript "air" when measured outside the glass. According to the notation adopted in the previous 
section no subscript indicates quantities measured in the glass, uppercase refers to the glass rest frame, lowercase to 
the frame comoving with the pulse at V ~ 0.66c as in fig. 2. 

We begin by studying the properties of the dispersion relation inside the glass, but outside the pulse (corresponding 
to the right panel of fig. 2). In fig. 6, loci of constant oj are plotted in the comoving frame (k x ,k y ). The implicit 
relation between k x and k y , in the comoving frame, at fixed uj, is directly provided by eq. (5). In the right panel 
the low momentum region is zoomed. The curves at fixed uj — w min and uj = w max are represented by blue and red, 
respectively, dashed lines. 

To explain the meaning and relevance of this plot, we first compare the section at k y = with the dispersion 
relation of fig. 2, right panel. For CL> m j n < uj < w max (colors from blue to red), there are four real solutions for k x , at 
fixed uj, represented by a single color between blue and red, and k y — (cf. the 4 solutions on the solid horizontal 
line of constant uj, in fig. 2). Then, when the frequency increases, the two rightmost solutions merge into a single one 
at uj = u; max (red dashed line) and disappear for uj > w max (colors from red to yellow). When, instead, the frequency 
goes down to uj — uj m - m (blue dashed line), the two central solutions collapse into a single point, around k x w 2/im~ 1 . 
Finally, for uj < w min (cyan curve), they completely disappear. In the full two-dimensional case, i.e., when k y ^ 0, the 
situation is more complicated and four solutions for k x may exist even for uj > w max , for sufficiently small but nonzero 
values of k y (see, for instance, the first red curve above the dashed red one). 

Noticeably, for all values of uj, there is always one solution on the leftmost branch. By comparison with fig. 2, right 
panel, this solution has negative rest frame frequency Q. The leftmost solution in fig. 2 lies indeed on the negative 
laboratory frequency branch, represented therein by a thin line. As noticed in the previous section, some quantum 
vacuum emission should therefore be present at any frequency uj. 

However, in order to be detected, phonons must be able to travel outside the glass reaching a detector at rest in 
the laboratory. To determine which photons can really reach the detector and what is their frequency and momentum 
in the laboratory, one must first make a Lorentz boost to the glass rest frame and then study how light exits the glass. 

Thus, in fig. 7 we report, in the rest frame of the glass {K x , K y ), the curves of fixed comoving frequency for the 
same values of uj as fig. 6, maintaining the same color correspondence. The relation between K x and K y , at fixed 
comoving frequency uj, is obtained by solving 



w = 7(V0[« - V K x ] 



(11) 
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Fig. 7. Loci of constant emission frequency uj (solid lines), equally spaced in log(aj), in the glass rest frame (K x , K y ), inside the 
glass. The dashed lines are to — o; m i n (blue), and to — tj max (red). Colors from cyan to blue refer to frequency ui < cj m i n , from 
blue to red to w m i n < uj < u; max , and from red to yellow to uj > cj ma x- In the right panel the low momentum region is zoomed. 
A mode does not undergo total reflection and can cross the interface between glass and air, if the corresponding wavevector 
(K x ,K y ) is contained in the green area. 



for fl and putting the result into eq. (2): 

c\Kl + Kl) = {uj/-! + vK x f nl{ujh + vK x ). (12) 

It is now possible to study whether a mode with a given uj, K x and K y can reach or not a detector placed outside 
the glass, by crossing the interface between glass and air. In fig. 7, the green area corresponds to rays that do not 
undergo total reflection and exit the parallelepiped through the lateral surface. Its boundary is the set of wave vectors 
(K x ,K y ) determined by imposing that K a \ ry vanishes in air, 

Ql r = c 2 Kl, x . (13) 

Since fl a i r = fl and K a \ rx = K x , from eq. (5) we obtain 

K 2 y = K 2 X [nl(cK x ) - 1] . (14) 

For the sake of completeness, we report in fig. 6 the corresponding locus of wavevectors describing photons that can 
cross the lateral surface of the parallelepiped in the comoving frame (k x , k y ). The procedure to determine this locus is 
analogous to the above one. However, the result depends on the sign of the rest frame frequency fl of the considered 
mode. The green area in fig. 7 is transformed into the green or blue area in fig. 6, when fl > or fl < 0, respectively. 
That is, if a mode has positive (negative) norm/laboratory frequency, then it does not undergo total reflection if the 
corresponding wave vector (k x , k y ) lies in the green (blue) area of fig. 6. 

It is now possible, either using figs. 6 or 7 to determine which modes can be observed outside the glass. Solutions 
on the leftmost and rightmost branches all lie in the white area, thus the corresponding modes cannot cross the 
interface parallel to the direction of the beam. Consequently, modes with u < oWm whose wave number lies on those 
branches (cyan curve), cannot be detected in the laboratory through the lateral surface 2 . We then consider modes with 
^min < uj < uj max (from blue to red curves). As shown in sect. 2.1, this is the frequency range where the scattering is 
described in the standard way by a 3 x 3 matrix. In this range, there are solutions with low momentum k, therefore 
quantum vacuum emission is expected to be stronger, by momentum conservation arguments. Unfortunately, only a 
tiny region of wave vectors with comoving frequency w m i n < uj < uj maK is contained in the green area. Finally, the 
largest part of the green area in the central panels is filled by modes with uj > cj max (from red to yellow curves) . Those 
photons, with very high comoving energy can cross the lateral surface and may be observed in the laboratory. 

After having determined which modes can exit the glass, we still have to understand how they look like in the 
laboratory, what is their wave number, and, in particular, what is their propagation direction. In fig. 8, we plot the 
loci of constant uj in air, in the glass rest frame (K a i TX , K a i Ty ), for the same values of ui as in figs. 6 and 7, with the 
same color legend. Since in the passage from glass to air, both the frequency and the x component of the wave vector 
are conserved, the dispersion relation in the comoving frame in air is 

c 2 (kl- 1IX + kl ily )=ujl il = uj 2 . (15) 

2 Of course, if the glass has defects, photons may be elastically scattered and exit the parallelepiped from the lateral surface, 
even if the x component of their momentum is initially large. 
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Fig. 8. Loci of constant emission frequency to (solid lines), equally spaced in log(oj), as observed in the laboratory reference 
frame where the glass is at rest (K^^, Kaiiy), outside the glass. The red dashed line is ui = w max . Colors from blue to red refer 
to u m in < uj < tJ max , and from red to yellow to uj > U) m&x . A mode can originate on the optical branch inside the glass if the 
corresponding wavevector (K a i TX ,K & i ly ) is contained in the green area. 

To determine the loci of constant uj in the laboratory reference frame in air, one must apply a Lorentz boost to this 
dispersion relation. However since this equation describes circles in the plane (k a i TX , k a i ry ) and a Lorentz transformation 
is affine, then loci of constant uj will appear as ellipses in the plane (K ailx , K aily ). Modes having wavevector in the 
green area can propagate inside the glass on the optical branch. The white area corresponds instead to modes that 
either do not belong to the optical branch (/2 a ir = fl < flgap) or are evanescent waves in glass. 

To conclude this section, we compare the loci of constant u> in the glass rest frame inside and outside the glass. 
The segments of curves lying in the green area of fig. 7 (rest frame in glass) are mapped onto the whole curves in fig. 8 
(rest frame in air). Instead, the segments of curves lying in the white area of fig. 7 do not have any counterpart in 
fig. 8 because of total reflection. That is, by looking through the interface parallel to the direction of propagation of 
the pulse, it is not possible to investigate the very interesting mode structure present inside the glass. 

Note that, in the laboratory frame, outside the glass, only a narrow region of wave vectors (K a i rx , K a i iy ) stays 
below the red dashed curve (uj — w max ), then corresponding to modes with comoving frequency smaller than uj max . In 
particular, no modes with uj < w max (from blue to red curves) propagate in air perpendicularly to the direction of the 
beam (if a ; r y = 0). Furthermore the blue dashed curve and the cyan curve of figs. 6 and 7 do not have any counterpart 
in fig. 7. This is not surprising because, as we already noticed, the blue dashed curve is not contained in the green 
area of figs. 6 and 7, implying that modes with comoving frequency uj < uj min cannot propagate in air, by crossing the 
lateral surface of the parallelepiped. We also noticed that the majority of modes that cross the lateral surface have 
frequency lo > w max (from red to yellow curves). A lot of them even propagate in the direction perpendicular to the 
beam, but unfortunately they cannot be photons produced in a standard analogue Hawking mechanism, because their 
comoving frequency uj is too high. 



4 Constraints from observations 

Here, we reverse the logic of the previous section, by determining the constraints put on the emission properties by a 
given observation in the laboratory. Let us assume that some observation is performed in the laboratory in a frequency 
range (f2 a i r i, ^ a ir 2 ) • What is the frequency window of emission u) compatible with those measures? The analysis is 
straightforward and the result is independent of the nature of the material. Since frequency is conserved when light 
crosses the interface ( f2 = i7 a ; r , uj — uj a i r ), radiation observed at frequency fl a i r in the laboratory must have been 
emitted at 

u} = 1 (V)[Q aiT -VK aixx ]. (16) 
Assuming again that the air refractive index is 1, one obtains 

uj = fl air j(V) [1-V cos # air /c] , (17) 

where 9 a i r is the angle of propagation in air with respect to the pulse direction. Letting cos# air vary from —1 to +1, 
the emission frequency window in the comoving frame is 



fa = flak i 7 00 [1 " V/C] ,UJ 2 = n ail2 j(V) [1 + V/C}) 



(18) 
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If measurements are performed by an instrument with frequency window (^ a iri, ^air2), only emission frequencies in 
the range (u^u^) can be explored. Using a pulse velocity of 0.7c, approximately typical of fused silica, uj\ = 0.4i7 a j rl 
and ll>2 — 2.4i7 a i r 2. In particular, photons produced in the comoving frame at low w cannot be detected. 

Conversely, let us assume that measurements are performed through a surface perpendicular to the direction of 
propagation of the pulse. In this case, radiation forming inside the glass a very small angle 9 with respect to the 
direction of the pulse impinges almost perpendicularly on the new interface. This radiation does not undergo total 
reflection and can be detected. The emission frequency uj in the comoving frame is then expressed in term of the glass 
rest frame frequency fl = J7 a i r as 



u = 1 {v)[n-VK x ] = 1 {v) nTVJn&n)n*/c 



K 2 

V 



V 



I f — \Jnl{Q &il ) - si 



sin 2 6» a 



(19) 



since the y component of the momentum is now conserved in crossing the interface. The upper (lower) sign applies to 
photons observed through the front (rear) interface. Since the group velocity of the pulse is not very different from 
the phase velocity at the pump frequency fl p , the pulse velocity is V ~ c/no(i7 p ). As a consequence, if f2 p is not far 
from the observation window (/2 a irij ^ a ir2)i n o(^air) ~ n o{^p) ~ C /V, and eq. (19) becomes 



1T\ 1 



V 2 



sin 2 alr 



(20) 



If measurements are performed through the front interface (0 < 9 < 90°, minus sign in the parenthesis), the emission 
window in the comoving frame compatible with observations in the range (/2 a irii Aiir2) is then (0,(7— l)^ a i r 2)- 
Interestingly, it goes down to very low frequencies, even if the observation window does not. This result implies that 
using a detector in a given frequency range (/? a iri> ^air2)j one can explore comoving emission frequencies in a very 
low frequency range. 

If measurements are performed through the rear interface (90° < 9 < 180°, plus sign in the parenthesis) the 
emission window compatible with observations in the same frequency range (i? a iri, ^ a ir2) hes instead at much higher 
frequencies ((7+ l)J7 air i, 27/? air2 ). 



5 Conclusions 

In the first part of the paper, we used simple kinematic arguments to identify under what conditions quantum vacuum 
radiation can be generated from a strong laser pulse propagating in a nonlinear Kerr medium such as fused silica. We 
identified three possible configurations: 

1) Horizon geometry. Even if the medium dispersion hinders the usual analogy with the propagation of a quantum 
field on a curved spacetime, frequency-dependent horizons can still be defined in a finite comoving frequency range 
w m in < lj < w max . Seen by a comoving observer, light within this frequency window can propagate only in one direction 
inside the pulse, while it propagates in both directions in the external region. In analogy with other analogue models 
based, e.g., on Bose-Einstein condensates, we expect that Hawking-like thermal radiation should appear for frequencies 
uj far enough from w m j n and w max . 

2) Horizonless geometry, fast pulse. If the pulse is fast enough, in the comoving frame modes can propagate only in 
one direction, both inside and outside the pulse. However, as negative-frequency modes exist also in this configuration, 
quantum vacuum emission is possible, but its spectrum is expected to strongly deviate from a Planckian law. 

3) Horizonless geometry, slow pulse. If the pulse is too slow, modes at frequencies oj below a threshold frequency oj t 
can propagate in the comoving frame in both directions both inside and outside the pulse, but a frequency-dependent 
horizon is present for oj t < uj < w max . Sufficiently far from oj t and w max , the emission is then expected to be Hawking- 
like as in 1). For all other frequencies, quantum radiation is still possible, yet with a different spectrum. 

On this basis, we are tempted to conclude that configuration 1) is the most suitable to observe the optical analogue 
of Hawking radiation. 

In the second part of the paper, we focus on configuration 1) to describe the propagation of quantum radiation 
from the pulse inside the nonlinear medium (where it is emitted) to the detector that is located in the surrounding 
air. In particular, we identified those modes that can cross the interface between glass and air. This analysis yields two 
important results: first, the Hawking-like emission within the w m i n < uj < w max window does not propagate outside 
the glass perpendicularly to the direction of propagation of the pulse. Second, the dispersive high momentum modes 
cannot cross a glass/air interface parallel to the pulse direction. This implies that the Hawking partners cannot be 
detected looking through this surface either. 

In conclusion, these kinematic arguments show that the optimal directions to observe Hawking-like emission in 
a realistic experimental setup as in [7] form a very narrow angle with the direction of the pulse. In addition to 
this, quantum radiation is anticipated to appear also in horizonless configurations. However, its spectral features are 
expected to strongly differ from those of standard analogue Hawking emission. 
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